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Running quark mass in two flavor QCD * 
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We present first results for the step scaling function op of the renormalization factor Zp of the pseudoscalar 
density. The simulations are performed within the framework of the Schrodinger functional with two flavors of 
O(a) improved Wilson fermions. The knowledge of <rp is required to compute the renormalization group invariant 
quark masses. We also study the performance of a variant of the HMC algorithm using two pseudofermion fields. 
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1. THE RGI MASS 

Lattice QCD is a theory which has as free pa- 
rameters the bare gauge coupling go and the bare 
current quark masses rrij , i = u,d,s,.... The 
hadronic scales like or tok are connected to 
the perturbative high energy regime of QCD via 
the running of renormalized couplings 



mi(/i) 



(1) 



where /i is the renormalization scale. At high 
energies the A parameter and the renormalization 
group invariant (RGI) quark masses Mi can be 
determined. 

In the following we assume that the renormal- 
ization conditions are independent of the quark 
masses themselves. The running of the renormal- 
ized couplings is described by the renormalization 
group equations (RGE) 



drrii 
dfi 



r(g)m l 



(2) 



The (3 and r functions are non-perturbatively de- 
fined if this is true for g and mj. Their perturba- 
tive expansions are 



-tibo + hg 2 + b 2 g 4 + ...}. 



(3) 



*Talk presented by F. Knechtli under the original title 
"Towards the computation of the strange quark mass with 
dynamical quarks". Work supported in part by the Eu- 
ropean Community's Human Potential Programme under 
contract HPRN-CT-2000-00145, Hadrons/Lattice QCD. 



r(g) ~ -9 2 {d + d{g 2 + ...}. 

9^0 



(4) 



RGI quantities P do not depend on the renormal- 
ization scale fj, 



li—P(n,g,mi) = 0. 



(5) 



Examples are the A parameter and the RGI quark 



x exp 



d.r 



1 1 bi 

(3{x) b x 3 b 2 x 



Mi = m l (26 5 2 )- do/2ho 



(7) 



x exp ■ 



dx 



T 0) 
P(x) 



dp 
b x 



where g — g(/i) and m; = rrii(p). These RGI 
parameters of QCD are defined beyond pertur- 
bation theory and their connections between dif- 
ferent renormalization schemes can be given in a 
simple and exact way. For these reasons A and 
Mi should be taken as the fundamental parame- 
ters of QCD. 

In this talk we present first results on the non- 
perturbative computation of the running of m(fi) 
in the Schrodinger Functional (SF) renormaliza- 
tion scheme with two dynamical flavors of O(a) 
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improved Wilson fermions. The strategy of our 
computation closely follows ref. We will 

give as a result an approximation to the flavor 
independent ratio M /m(p) as a function of /i/A. 

2. THE RUNNING MASS 

In continuum QCD a renormalized mass in the 
SF scheme is defined through the PCAC rela- 
tion which involves the renormalized axial cur- 
rent (Ap/j^x) and the renormalized pseudoscalar 
density Pr(x) 



Pr(x) 



(ffii + rfij)PR, 
Zptpi(x)j 5 ^j(x) . 



(8) 
(9) 
(10) 



The renormalization constant Za is fixed by 
Ward identities. The renormalization constant 
Zp is defined as 



Z P (L) 



at rrii 



0. 



(11) 



ML/2) 

For notation and details of the calculation we re- 
fer to [0J|]. The correlation functions /1 and fp 
are computed at zero bare current quark masses 
rrn. The renormalization scale in the SF scheme 
is identified with the inverse spatial box size /i = 
1/L. Keeping the physical size L constant can 
be realized by keeping the renormalized coupling 
g(L) fixed. 

A rigorous definition of the renormalized mass 
can be given in the lattice regularization of QCD 
by 

Za(9q) 



m,i(fj,) = lim 



: m i(go) 



(12) 



a— >o Zp(ga, L/a) 

The running of the mass can be described by the 
flavor independent step scaling function up(u) 

m(n) = Z P (g ,2L/a) 



m(/i/2) a->o Zp(g ,L/a) 



a P (u). (13) 



u=g 2 (L) 



Together with the step scaling function for the 
gauge coupling a(u) = g 2 (2L) we can solve a 
system of coupled recursions: k = 0, 1, 2, ... 

MO-9 2 (imax)= 3.3000 

u k = g 2 {2- k L max ) 



(14) 
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Figure 1. The running of the renormalized quark 
mass. 



v = 1 

m(/x m i„) 



(15) 



t'k 



m(2> min ) 



2L n 



The largest coupling uq defines the reference scale 
L mM . At high energies contact with perturbation 
theory can safely be made and the RGI parame- 
ters are determined from eq. @ and eq. (Q) using 
the perturbative 2-loop t function and 3-loop (3 
function. 

We computed the lattice step scaling function 
Sp(u, a/L) for 6 couplings in the range u = 
3.33 ... 0.98 on L/a = 6 and L/a = 8 (requir- 
ing also 2L/a) lattices. We interpolate Ep by the 
Ansatz 

E P (u,a/L) = 1 - ln(2)d u + p 2 u 2 + p 3 u 3 (16) 

with fitted parameters p 2 and p%. Using it in the 
recursion eq. (|l5| ) we obtain the following approx- 
imations for 



M 



m(/i min ) 



1.202(13) 
1.236(15) 



L/a = 6 
L/a = 8 



(17) 



(contact with perturbation theory is made at 
the scale 2 _8 L max ). Combining with recursion 
eq. (|T^), from which we get in the continuum limit 
ln(AL max ) = -1.85(13) @, we plot in Fig. the 
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Figure 2. Performance of the HMC algorithm us- 
ing one or two psedofermion fields (PFs). 
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(18) 



as a function of /i^ = 2 /i m i n , together with the 
pertubative curve using the 2-loop t function 
and the 3-loop j3 function. The errors of the 
points come from the statistical errors in the coef- 
ficients Vk ■ The overall uncertainties in the quan- 
tities AL max and M /m(fi m i n ) are not shown. The 
points show a dependence on the lattice spacing 
but for an extrapolation to the continuum limit 
we need additional simulations on finer lattices. 



3. ALGORITHM STUDY 

In our simulations we compared the computa- 
tional cost of the HMC algorithm for the standard 
action and for an action which uses two pseud- 
ofcrmion fields (PFs) g. We denote by Q the 
even-odd preconditioned hermitian Dirac opera- 
tor for 0(a) improved Wilson fermions. The ef- 
fect of det(Q 2 ) is taken into account by two PFs 
4>\ and (f>2 with actions 



Sf 2 



(Q 2 



4 p- 2 



p ) 

Q- 2 



(19) 
(20) 



This splitting of the action is valid for any real 
parameter p. We take 

P = fA mi n(0 2 )A max (g 2 )) 1/4 , (21) 



which minimizes the sum of the condition num- 
bers in eq. ( p^ ) and eq. (|20|). This is expected to 
reduce the fermionic force allowing for larger step 
sizes in molecular dynamics. 

The implementation of the Hybrid Monte Carlo 
(HMC) algorithm to simulate the system de- 
scribed by eq. ( |i~9| ) and eq. ( pfj| ) is straightfor- 
ward. Simulations with the same physical pa- 
rameters show that the step size of molecular 
dynamics with two PFs can be doubled with re- 
spect to the standard case with one PF by keep- 
ing the same acceptance (optimal around 80%). 
In Fig. H we show the performance of the HMC 
algorithm on the APEmille machine, 1 crate of 
the machine consists of 128 nodes and sustains 68 
GFlops peak. The computational cost to achieve 
1% error in Zp at constant value of ~g 2 k 2.48 is 
almost a factor 2 lower for two PFs than for one 
PF on lattices L/a = 12 and L/a = 16. This gain 
in performance is important in view of the next 
simulations on lattices L/a = 12 and 2L/a = 24 
needed to obtain ap(u) in the continuum limit. 
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